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Abstract 

We consider an interaction representation in the Boltzmann field theory. It de- 
scribes the master field for a subclass of planar diagrams in matrix models, so called 
half-planar diagrams. This interaction representation was found in the previous pa- 
per by Accardi, Volovich and one of us (I. A.) and it has an unusual property that one 
deals with a rational function of the interaction Lagrangian instead of the ordinary 
exponential function. Here we study the interaction representation in more details 
and show that under natural assumptions this representation is in fact unique. We 
demonstrate that corresponding Schwinger-Dyson equations lead to a closed set of 
integral equations for two- and four-point correlation functions. Renormalization of 
the model is performed and renormalization group equations are obtained. Some 
model examples with discrete number of degrees of freedom are solved numerically. 
The solution for one degree of freedom is compared with the planar approximation 
for one matrix model. For large variety of coupling constant it reproduces the planar 
approximation with good accuracy. 
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1 Introduction 



There is an old problem in quantum field theory of finding the leading asymptotics in 
N X N matrix models for large in realistic space-time. Its solution may have important 
applications to hadron dynamics since the large N limit in QCD {N is the number 

of colors) enables us to understand some phenomenological features of strong interaction. 

In the early 80-s it was suggested |0 that there exists a master field such that 
correlation functions of this field are equal to the large N limit of invariant correlation 
functions of a matrix field M{x), 

The problem of the construction of the master field has been discussed in many works, 
see for example 0-[|18|. Gopakumar and Gross [|19] and Douglas have constructed 



the master field for an arbitrary matrix model in terms of correlation functions by using 
methods of non- commutative probability theory |^-|2^. There has been a problem of the 



construction an operator realization for the master field without knowledge of correlation 
functions. Recently this problem has been solved in and it was shown that the master 
fields satisfy to standard equations of relativistic field theory but fields are quantized 
according to a new rule. These fields have a realization in the free (Boltzmannian) Fock 
space. 

Quantum field theory in the Boltzmannian Fock space has been considered in [EB, IH, 



pq , |2^. Some special form of this theory realizes the master field for a subset 

of planar diagrams, for the so called half-planar diagrams. This construction deals with 
the master field in a modified interaction representation in the free Fock space. This new 
interaction representation involves not the ordinary exponential function of the interaction 
but a rational function of the interaction and correlation functions are given by the formula 

< $(xj...<i>(xi) >= (o|0(xJ...0(xl)|^]), \n) = n\o) 

^= ,^0 f.y (1-2) 
where (l){x) is a field in the free (Boltzmannian) Fock space, 

(f){x) = (f)-{x) + (f)-^{x), (1.3) 

satisfying the relation 

r(a:)0+(y) = A(a;-y), (1.4) 

where A{x — y) is the free propagator and 0~|O >= 0. 

In this paper we study the Boltzmann field theory with correlation functions ( p..2|) and 
we show that the Boltzmann field theory gives an analytical summation of half-planar 
diagrams. 

We prove in pure algebraic way that the correlation functions ( [L.2D satisfy a closed set 
of Schwinger-Dyson-like equations. Moreover, starting with an arbitrary function of the 
interaction Vt = Q{Sint{(p)) and using natural assumptions we show that the form ( |1.2|) of 
^{Sint) is in fact an unique one which admits Schwinger-Dyson-like equations. We will 
call these equations the Boltzmannian Schwinger- Dyson equations. 
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We will see that there are an essential simplifications in the Boltzmannian Schwinger- 
Dyson equations. Namely, among the Boltzmannian Schwinger-Dyson equations for n- 
point Green's functions for n greater then the degree of Sint there are equations which 
relate n-point Green's functions only with A;-point Green's functions for k < n. This is 
a distinguish feature of the Schwinger-Dyson equations for the Boltzmann fields and it 
drastically simplifies the situation. In particular in the case of quartic interaction one 
has a closed set of equations for two and four point correlation functions. Let us stress 
that this system of equations is exact one but not a truncation of the original system 
and it does not assume any approximation for correlation functions ( |1.2|) . The origin of 
this property of the Boltzmannian Schwinger-Dyson equations is related with a specific 



features of the Wick theorem for the free algebra (|1.3|) . Note that in the operator form 
the Boltzmannian Schwinger-Dyson equations contain projectors on the vacuum state in 
the interaction terms. 

The Boltzmannian system ( |1.2|) can be considered as a non-trivial approximation 
to the planar correlation functions. Note in this context that in all previous attempts 
of approximated treatment of the planar theory some non-perturbative approximations 
were used |T^, |TB|, ^ . Topologically diagrams representing the perturbative series of the 
Boltzmann correlation functions ( |1.2| ) look as rainbow (or half-planar) graphs of the usual 
diagram technique for matrix models |^ . This is a raison to call the Boltzmann correla- 
tion functions the half-planar correlation functions and the corresponding equations the 
half-planar Schwinger-Dyson equations. Comparing three sets of the Schwinger-Dyson 
equations, full, planar and half-planar ones, one can say that the operator form of the 
full Schwinger-Dyson equations is the simplest one. Indeed, the planar Schwinger-Dyson 
equations contain projectors in the Schwinger terms and the half-planar Schwinger-Dyson 
equations also contain projectors in the interactions terms. However just due to the pres- 
ence of additional projectors it turns out that the half-planar Schwinger-Dyson equations 
admit analytical investigations. 

We will find explicitly the correlation function ( p.. 21 ) for the quartic interaction in 
the case of one degree of freedom. We compare numerically the two-, four- and six- 
point half-planar correlation functions with the corresponding planar correlation functions 
for the one matrix model. For large variety of the coupling constant the half-planar 
approximation reproduces the planar approximation with good accuracy. However the 
half-planar approximation does not reproduce asymptotics of correlation functions in the 
strong coupling regime. To reproduce the strong coupling regime one has to consider 
more large class of graphs 



For the case of finite number of degrees of freedom and a special action with a diag- 
onal quadratic term it is also possible to get an analytical solutions of the Boltzmannian 
Schwinger-Dyson equations. For more complicated case of a non-diagonal quadratic ac- 
tion and quartic interaction the system of equations for two- and four-point correlations 
functions is reduced to an algebraic system of equations that admits a numerical solu- 
tion. This system describes the half-planar approximation of the two-matrix model which 
has been analyzed recently by Douglas and Li |^ using the language of functions of 
non-commutating variables. 

An analytical investigation of the correlation functions ( |1.2| ) is also possible for D-di- 
mensional space-time. In the case of quartic interaction one has a closed set of integral 
equations for two- and four-point correlation functions. For four-point correlation func- 
tions we get a Bethe-Salpeter-like equation. A special approximation reduces this system 
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of integral equations to a linear integral equation which has been considered in the 
rainbow approximation in the usual field theory. 

For four-dimensional space-time the perturbative expansion of correlation functions 
( |1.2| ) has ultraviolet divergences. To remove these divergencies we apply R-operation 



33| . We show that an application of R-operation is equivalent to an introduction of 
counterterms to the interaction Lagrangian. As usual this fact permits to write down the 
renormalization group equations. Because of special features of the Boltzmann field theory 
the renormalization group equations are different from the ordinary renormalization group 
equations. There is also a difference in the explicit formula for the beta-function in terms 
of renormalization constants. 

The paper is organized as follows. In Section 2 we deal with one degree of freedom, in 
Section 3 we consider finite number of degrees of freedom and the last section is devoted 
to continuous fields in D dimensions. We start Section 2 by reminding specific feature of 
free Fock space for one degree of freedom. In subsection 2.2 we show in pure algebraic way 
that the correlation functions ( p.. 21 ) satisfy a closed set of Schwinger-Dyson-like equations. 
Then we study a question of uniqness of the form ( |1.2|) for Vt. In 2.3 we present an operator 
form of the Boltzmannian Schwinger- Dyson equations. In 2.4 we solve explicitly these 
equations. In subsections 2.5, 2.6 and 2.7 we present a generating functional for correlation 
functions for different forms of interaction. In subsection 2.8 we compare the half-planar 
approximation with the exact answer for the sum of the planar diagrams for the one matrix 
model. In Section 3 we consider the same questions for finite number of degrees of freedom. 
The goal of Section 4 is to demonstrate that the infinite set of the Schwinger-Dyson 
equations for Boltzmann correlation functions (|1.2|) in D-dimensions allows a reduction to 
a finite system of integral equations for lower correlation functions. In subsection 4.2 we 
consider the quartic interaction and present a system of integral equations for two- and 
four-point correlation functions. In 4.3 we study the renormalization of the correlation 
functions ( [L.2|) and derive the renormalization group equation for D = A. In Appendix a 
calculation of combinatoric factors for planar correlation functions is presented. 



2 One Degree of Freedom 

2.1 Boltzmannian Fock Space and Free Correlation Functions 

In this section we consider the free correlation functions in the Boltzmannian Fock space 
in zero-dimensional space-time. Let us start from the general definition of Boltzmannian 
Fock space, then we will restrict ourself to zero-dimensional case. 

The free (or Boltzmannian) Fock space Ti over the Hilbert space H is just the tensor 
algebra 

Creation and annihilation operators are defined as 

a^{f)fl ® - ® /n = / ® /l ® - ® fn 

® - ® fn =< /, /l > ®/2 ® - ® fn (2.1) 

where < f,g > is the inner product in H. We have 

aif)a+ig)=<f,g>. (2.2) 
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We shall consider the simplest case H = C. One has the vacuum vector |0), 

a|0) = (2.3) 

and n-particle states 

|n) = (a+)"|0). (2.4) 
The master field in zero-dimensional case is 

= a + a+, (2.5) 

where a and a"*" satisfy the following relation 

aa+ = 1. (2.6) 

This algebra has a realization in the free (or Boltzmannian) Fock space. 

One can formulate the Boltzmannian Wick theorem which is an analog of the ordinary 
Wick theorem. For this purpose let us define the normal product of operators ( |2.5| ) in the 
following way 

^ (a+)'"a'^-'". (2.7) 

m=0 

Note that ( p.7| ) is different from the normal product of n-t\i power of ordinary field 

if = a + a'^ 

in-m)\m\ 

where [a, a^] = 1. For an arbitrary product of creation and annihilation operators we 
define the normal product as 

: {a+y''{a)"''{a+y'^{af'\..{a+p{a)"" : = 

j:{a^r''{af'=6^^^,6M,,o, (2.8) 

k=l 

k ^ I k-1 ^ I 

1=1 i=k i=l i=k+l 



where 



and 

' 1, n = 0, 
0, n ^ 0. 

Let us call a contraction of two operators an admissible one if it relates neighboring 
operators. The Boltzmannian Wick theorem states that the product of operators ( p.8| ) 
is equal to the normal product of these operators plus sum of normal products of these 
operators with all admissible contractions. In particular, 

min{n,m) 

: 0" :: 0"* :=: 0"+™ -+12 ■ 0''+"'"^'' : • (2.9) 

fe=i 

A free n-point Green's function is defined as the vacuum expectation of n-th power of 
master field ( p3|) 

Gi") = (O|0"|O). (2.10) 
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As it is well-known, the Green's function ( ^.1U| ) is given by a n-th moment of Wigner's 
distribution 0, ^ 

2n J-2 



:2.ii] 



and 



2n 



Cr,, where the Catalane numbers 



(2n)! 



n!(n + 1)! 



(2.12) 



In an algebraic origin of the representation ( p.ll|) has been found. It is related with 
the fact that the algebra ( |2.6| ) is isomorphic to the quantum semigroup SUq{2) with g = 



and ( |2.11|) is the Haar measure on SUq{2). 

The representation ( |2.11 ) can be also obtained as a solution of the Schwinger- Dyson 
equations. To get these equations from definitions ( |2.6|) - (|2rT0|) let us calculate (|2.10|) using 
the Boltzmannian Wick theorem. We present the result of an application of the Wick 
theorem graphically (Fig. 1). We draw the n-th power of as n points lying on a line and 
the contractions as upper half-ovals. According to the relation aa^ = 1 all points may be 
connected only by non-overlapping upper half-ovals. Considering all possible contractions 
of the first point from the left and taking into account that all odd correlation functions 
are equal to zero we have 

n 

r^(o) _ r^o) MO) ^9 , n\ 

^2n — ^2m-2^2n-2m- l^-J-Jj 

m=l 

Equation (|2.13|) leads to an algebraic equation ^6] for a generation functional 



n=0 



Indeed, multiplying (|2.13|) by j^" and making the summation over n we have 



(2.14) 



Eryi0)-2n _ ^(0) ^(0) .2n 

^2nJ ~ 2^ 2^ '^2m-2^2n-2mJ ■ 



n=l 



n=l m=l 



Using 

oo n oo oo 

n=l m=l m=l n=m 

and making the change m — 1 m, n — m —>■ n, we get the equation 

oo oo oo 

E^(0) -271 _ Y^ Y^ ^(0) ^(0) .2n+2m+2 

^2nJ ~ 2^ 2^ ^2m^2nJ i 
n=l m=0 71=0 

which in terms of generation functional can be written in the form 
Under the initial condition Z'^'^)(0) = 1 the solution of equation (|2.15|) is 



3] 



1- y/T^lp 

2p 



- ^2n)^ 
t,n\{n + l)\^ 



(2.15) 



:2.16) 
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One has the integral representation for Z'^'^^j) 

ITT J-2 1 — A? 



[2.17) 



that is in agreement with ( |2.11 



Note that in the above formulae j may be considered as an operator. If one assumes 
that j satisfies a non-standard differentiation rule 



d 



3' 



then one can deduce the following differential equation for Z{j) [0: 
To reproduce ( p. 161) one has to assume the following initial conditions 



zW(o) = i, (4^(°))(o) = o. 

dj 



[2.18) 



[2.19) 



(2.20) 



In what follows we shall also use the following modified form of equations (|2.13|) 



fc-i 



r^(0) _ ^(0) ^(0) ^(0) ^(0) /r, 9. ^ 

'~^2n — "-^2A;-2m-2'-^2n+2m-2fc '-^2m-2fc'-^2n.+2fc-2m-2 l^'^-^y 



m=0 



m=k 



for any k = 1, 2, n. 
index. Indeed, making the change m —>■ m' = k — m in the first sum of ( p.21|) , we get 



Equations ( |2.21| ) are obtained by a simple change of summation 



k 

/ . '-^2m'-2"-^: 



2n-2m'- 



m'=l 



Introducing the new summation index in the second term in right hand side of (|2.21|) as 



n + — m we get 



n'G. 



(0) 



'-^2n-2m"-^2m'-2' 
m'=fc+l 

SO, (|2.21|) is equivalent to ( p.l3|) . Equations (|2.21| ) are the Schwinger-Dyson equations 
written for the case when one traces down for contractions of 2fc-th operator in the vacuum 
expectation (00... 0). This is illustrated on Fig. 2. The same equations are obtained when 
one traces down for the contraction of 2k + 1-th operator (Fig. 3). 

Equations (|2.21|) are 0-dimensional analogues of the D-dimensional Schwinger-Dyson 
equations written for the case when the operators (— A^^^fe + "^^) ("^2^2^+1 + '^^) 
applied to a Green's function G2'n{xi, ...,Xn)- 

2.2 Boltzmannian Schwinger-Dyson Equations 

In the ordinary Euclidean formulation of Bose or Fermi 0-dimensional quantum field 
theory the interacting Green's functions have the form 



^(Euci.) _ (o|(^"e-^^"*(^)|0). 



(2.22) 
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A natural generalization of this formula to the case of Boltzmann statistic is 

g., = {Q\ct>^n{sU4>)m, (2.23) 

where VL{Sint{<P)) is not necessary an exponent but some analytical function of the La- 
grangian. First of all we will consider the function VL{Sint{4>)) of special form |^ 



^{S^nt) = (2.24) 

and we will derive the Schwinger-Dyson-like equation for interacting Green's functions 
(|2.23|) . Then under some natural assumptions discussed below we will prove that the 
form ( p^.24| ) is unique. 

We will consider the case of quartic interaction Sint = g4>'^- The generalization for an 
arbitrary polynomial interaction is straightforward. 

Using the Taylor expansion for Q{g(f)'^) 

oo 

nig<p') = y: g""^^^'"", (2.25) 

m=0 

where Qm = (—1)™ for ( |2.24| ), we write the Green's functions Qn in terms of the free 
Green's functions Q^^: 



Qn=Y: ^7"^fi^(O|0"+^'"|O) = r^mGl^Um- (2-26) 

m=0 m=0 

Using the free Schwinger-Dyson equations in the form 

G'^' = E e^-i^ri-.-i + E k = 1, ...,n, (2.27) 

1=1 l=k+l 

(that coincides with (|2.21| ) after setting Qn = for odd n) we get 

oo k—l n+im 

= E 5'™'^m[E ^fc-i-l^«+n+4m-fc-l + E ^ l-k-l^ n+Am+k-l-l\- (2.28) 
m=0 1=1 l=k+l 

Writing the sum in the second term of (|2.28|) as 



n+4m n n+Am 

E = + T. (2-29) 

l=k+l l=k+l l=n+l 

and making the change of the summation index I — n —>■ I in the second sum of ( p. 29] ), 
one obtains 

k~l n 

= E Gk-l-lGl+n~k-l + E Gl~k-lGn+k-l~l + 
1=1 l=k+l 

oo Am 

E a'^^m E Oll'li-k-Ak-i-i = Ji + J 11 + J in- (2.30) 

m=0 1=1 
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Using the formulas 

4m m—l 



E /(O = E [/(4P + 1) + /(4p + 2) + /(4p + 3) + /(4p + 4)], 

oo m—l oo oo 

E E/Kp) = E E /Kp) (2-31^ 

m=0 p=0 p=0 m=p+l 



and making the change m ^ m + p + 1 in the second hne of ( p.31|) one can rewrite Jm 
in the following way 

oo oo 

■^-f-f-f ~ 2^ y ^'m.+P+lL^n-/c+4p^4m+A;+2 "T ^n-fc+4p+1^4m+/c+l"r 

p=0 m=0 

^(0) ^(0) ^(0) ^(0) n .2 321 

^ ='n-fc+4p+2»4m+fc »n-fc+4p+3^4m+A;-lJ ) l^Z.OZJ 

so equation ( |2.3CI| ) takes the form 

k-l 



- E Sk^-l-l^l+n-k-l + E Sl-k-lSn+k-l-r 
1=1 l=k+l 



y ^^rn+p+l[yn-k+4p^4:m+k+2 ' ^n-k+ip+l^im+k+l' 

p=0 m=0 

C(0) (-(0) ,M0) MO) 1 /2 33) 

='n-fc+4p+2^4m+fc ' ^n-fc+4p+3^4m+fc-lJ • I^^.OO; 

Note that two first sums in ( p.33| ) are nothing but Schwinger's terms. Since the 
coefficients Qm satisfy the relations 

^n+m ^n^m; (2.34) 

we have a possibility to represent the third sum as a combination of interacting Green's 
functions ( p.23 ). With Q in the form (|2.24 ) equations ( 2.33 ) take the form 

Gn = —y[Qn~kQk+2 + Qn-k+lGk+1 + Gn-k+2Qk + ^n-fc+3^+fc-l] + 
k—l n 

E ^fe-i-l^i+n-fc-l + E Gl-k-l^n+k-l-l- (2.35) 
1=1 l=k+l 

Equations ( |2.35| ) are the Schwinger-Dyson-like equations for zero-dimensional Boltzmann 
theory with quartic interaction. We will also call these equations the Boltzmannian 
Schwinger-Dyson equations. 

We see that the terms describing an interaction are quadratic on the correlations 
functions. The sum of indices in all these terms is constant and depends on the power 
of the interaction. One can prove that under the assumption of such structure for the 
interaction terms in the Schwinger-Dyson-like equations the form of fl is unique. In fact, 
one assumes that the form of Q does not depend on the special choice of the interaction 
and Q{0) = 1. These assumptions mean that in particular for the linear interaction 
Sint = y(p one has the following Schwinger-Dyson-like equations 

Qn = —yGn-iGo + Schwinger's terms, n > 0. (2.36) 
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Let us derive that the form ( |2.24| ) follows from ( ^.361 ). Expanding the both sides of ( ^^.361 ) 
on power coupling we have 

oo oo oo 

g'^^mQ^nlm = '9 ■ E d'^kG^'^ + Schwiuger's terms, (2.37) 



m=0 



m=0 



k=0 



Using equation ( p.27|) and comparing the coefficients in front of the same powers of the 
coupling constant we obtain 



(2.38) 



k=0 



For m = from ( p.38|) we have fl^ + fli = and taking into account fio = 1 we get 
Qi = —1. Assuming 

= (-1)"^ (2.39) 

for m < tuq and using equation (p.38|) one can prove that this relation is also hold for 
m = mo, that gives by induction the prove of ( p.39| ). 

From the above consideration we see that to get a closed set of Schwinger-Dyson 
equations in the free Fock space for an interacting case one has to deal with a non- 
Heisenberg dynamics and the Green's function in the interaction representation are not 
represented as exponential function of the interaction but are given by formula p.24| ). 



2.3 An Operator Form of the Boltzmannian Schwinger-Dyson 
Equations 



Let us write down the week operator form of equation (|2.35|) . By the week operator 
form we mean the relation that is hold after averaging with the function ( |2.24| ). Let first 
consider the case k = n. Then equations (|2.35| ) take the form 



(2.40) 



The left hand side of ( |2.40| ) represents the term 00" ^. The last term in the right hand 
side of ( |2.40| ) is the Schwinger term and it can be represented by the following commutator 



Z TT, 



where vr is an element of the following algebra 

= ilO >< 01. 



Indeed, after performing the average of (|2.41|) with Vt we get 



jt-2 



- I < 0|[7r, 0"-l]^]|O >= '3f'3n-l-2 



(2.41) 



(2.42) 



^2.43) 



l=Q 



The sum in the square brackets in ( p.40| ) can be represented as an average of the product 
of 0"~^ with the following sum of operators 



\VL >< 0\+(h^\Q >< 016 + SQ >< 016^ + \Q >< 01 



(2.44) 
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Note that these terms can be represented as 



5<p 



502 



...+ 



int 



|fi)(o| 



if one assumes that non-commutative derivatives 



in-l 



(2.45) 



(2.46) 



are imphed in (p.45|) . Therefore, equations (|2.40|) is the equation for the average 

< O|^0"-l|^] >= 0, (2.47) 
where the operator £ can be represented as 

^ = ^ + V.., + 5 (2.48) 
o(p[x) 

with 

5= [tt, ■ ]. (2.49) 
For k < n equations ( p.35|) can be also represented as the overage 

< 0\(j)''S(j)'''''-^Q >= 0, (2.50) 

where one has to use another form for S. 

2.4 Solution of the Boltzmannian Schwinger-Dyson Equations 

The distinguish feature of equations ( p.35| ) is that for n > A and 2 < k < n — 1 the right 
hand side of (p.35|) does not contain the Green's functions Qm with m > n. This fact 
permit us to write down a closed system of equations for any Qn, n > A. Let us study 
this questions in more details. 

Equations ( |2.35| ) for /c = 1 and k = 2 take the form 



n 

Qln = 5Z ^2i-2^2n-2i " 5'^2n^2 " 5'^2n+2^0) 
1=1 

n-1 

Q2n = X! ^2«-2^2n-2i + ^2n-2 " fl'^2n-2^4 — fl'^2n^2, 
1=1 

where we use the fact that Q2n+i = 0. Setting n = 1 in ( p.51|) and n = 2 in ( p.52| 
obtain the system of two equations for three unknown functions Go{g), G2{g) and ^4 



we 



equations for three unknown functions Go{g), G2{g) and Q. 
G2 = Go — fi'^2^2 — gGiGo, 

Gi = 2G2 - 2gG2GA- (2.53) 
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To have the third equation let us note that there is an additional set of equations for 
Green's functions 



G2n — ~ 5'^2n+4; 



that follow from the identity 



(o|- 



k2n 



:|O) = (O|0^"- 



9<P' 



,2n+4 



-|0). 



For the vacuum Green's function Qq equation ( |2.54| ) has the form 

^0 = 1- gGi- 

The system of equations ( p.53| ) and ( |2.56| ) lead to the following equation for Qq 

AgQt + - 1 = 0. 
Taking into account the condition Qq = 1 ioi g = Q have 



(2.54) 



(2.55) 



(2.56) 



(2.57) 



-1 + yi + IQg 

8g 



:2.58) 



The Green's functions Q2 and Q4 can be expressed in terms of Qq: 



2ggo 
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(2.59) 



High-point Green's functions can be found from (|2.54| ) after taking into account ( p.58|) 
and ( P3D|) . 

2.5 Generating Functional 

A generating functional for interacting Green's functions is 



^(j,^?) = E^2„(^7)f". 



(2.60) 



n=0 



From ( |2.11| ) one has the following integral representation for Z{j,g) 



[2.61) 



and the integral ( |2.61| ) can be calculated explicitly. However it is instructive to get 
an explicit derivation of Z{j,g) in a pure algebraic way, since just this method can be 
generalized to the D-dimensional space-time. 
From (|23l| ) it follows 



f[ZU) - Go] = fZ^'\j)ZU) - gfg^lZU) - Go] - g\m - G2f - Go]Go (2.62) 



and 



fg,{i + 2gg2) + ggi 
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2go{f + ggi) 



Using ( p.54|) one can write the equation for Z{j) in an ahernative form 



'z{j)=j'z^'\j)-g[z{j)-g,f-go\ 



(2.63) 



(2.64) 



and therefore 



fZ^){j)+gg2f + gg, 
3^ + 9 



2ggoZ^'\j) + gf - ggof + 2g'gi 



2ggo{3^ + g) 

Taking into account (|2.15|) and ( p.57|) one can verify that (|2.65| ) coincides with ( |2.63| ). 



2.65) 



2.6 Correlation Functions without Vacuum Insertions 

To calculate the Green's functions without vacuum insertions G„ note that they satisfy the 
same Schwinger-Dyson equations, but one has to use another "initial" condition, namely 
Go = 1- The system of equations for 2-, 4-, and 6-point Green's functions is 

G2 = 1 — gG2G2 — gG^, 
G4 = 2G2 — 2gG2G4j 



Gq — G2 + 2G4 — gG\ — (7G2G6. 
From the first two lines of ( |2.66| ) one gets the following equation for G2: 

2g'^Gl + ^Gl + G2 - 1 = 0. 

Introducing a new variable G2; 



Go 



Go 



2g 



we have 



Gl-G2-4:g = 0. 



(2.66) 
(2.67) 

(2.68) 
(2.69) 



The solution of equations (|2.69|) has a different form in dependence on values of g. 
Using the Cordano formula we have 



G2 = -^Ree^t( Jl - 108^2 _ iQ,/sg)l 
v3 



(2.70) 



for g < and 



G2 = -^{(^^9 + /lOsT^)^ + {QV3g - 7108(72 -1)^; 



(2.71) 



for (7 > -tI . Therefore 
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G2 = ^[-^Ree*t(v'l - 108^2 _ ,Q^gf^ _ 1] 
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1-3^ + 16/- 105^^ + 768^^.. (2.72) 

for g < ^ and 

G2 = ^^[(6v^^? + \/l08(72-l)i + {QV3g - ^W-l)! - v^] (2.73) 
for g>4. 

An equation for a generating functional for Green's functions without vacuum inser- 
tions Z{j) follows from (p.62|) by setting Gq = 1: 

f[Zij) - 1] = fZ^^\j)Zij) - gfC^iZij) - 1] - g[Zij) - G^f - 1]. (2.74) 



From (g:74D we find 



g+ngG2 + i)-fzi%jy ^'-''^ 



where Z^^\j) and G2 are given by ( p.l6|) and ( |2.72|) , (|2.73| ), respectively. 



The generating functionals Z{j,g) and Z{j,g) are related by the formula 

Z{j,g) = Goig)Z{j,gGo{g)), (2.76) 
which follows from the relation 

Gnig)=Goig)GnigGoig)). (2.77) 

The last relation is true due to the Boltzmannian Wick theorem. It is easy to check that 
Z{j,g) and Z{j,g) given by (|2.63|) and ( p.75| ) satisfy the relation ( p.76| ). Indeed, from 



(|Z75|) it follows that 

gGl+fil + 2gGoG2{gGo)) 



Goig)Z{j,gGoig)) 



gGo + f{l + gGoG2{gGo) - fZi^){j) ' 

Replacing in the right hand side of this relation GoG2{gGo) by G2 and taking into account 
(|2.62|) we get the expression for Z{j,g). 



2.7 Correlation Functions for Normal Ordered Interaction 

In this subsection we will consider the Boltzmann theory with a normal ordered interaction 
: 0^ :. Let us denotes by Gn the n-point Green's function without vacuum insertions in 
such theory. According to the definition of Boltzmannian normal product we have 

: (j)^ ■= 04 _ _ 2 (2.78) 

and so 

G2n = (0102" —r^ )|0). (2.79) 

^ 1 + ^(04-302-2^' ' ^ ' 

Let us write down the Schwinger-Dyson equations for Green's functions G„. Repeating 
all steps as for derivation of ( |2.51| ) and (|2.52|) we get 



n 



G2n — X/ ^2Z-2^2n-2« — gG2nG2 — gG2n+2Go + 3gG2n, (2.80) 



1=1 
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n-1 



G2n — ^ G2i_2G2n-2l + G2n~2 — gG2n-2Gii — gG2nG2 + 3gG2n-2G2- (2-81) 



1=1 



Putting n = 1 in ( p^.80D , n = 2 and n = 3 in ( |^.81| ) we obtain the following equations 
for (^2, and Gq 

G2 = l-gG2G2-gG, + 3gG2, 
(j4 = 2G2 — 2gG2G ^ -\- 3gG2G2- 
Gq = G2 + 2G4 - gGl - gG2Ge + 3gG2G^. (2.82) 
From the first two lines of (|2.82|) one can derive the equation for G2 



2g^Gl + 3g{l - g)Gl + (1 - 3g)G2 -1=0. (2.83) 
We restrict ourself by writing the perturbative solution of ( p.83| ) up to order g'^: 

G2 = I + /-3/ + 9/ + ... . (2.84) 

2.8 The Half-Planar Approximation and the Two-Field Boltz- 
mann Theory 

In the previous subsections we have investigated the interacting Boltzmann theory itself. 
In this subsection we are going to compare the Boltzmann theory and the planar ap- 
proximation for the one-matrix model. Green's functions for the one-matrix model in the 
planar approximation are defined as 

n2n{g) = Ji^^]^^ / t^Mtr(M2")exp[-^tr(M2) -A-tr{M% (2.85) 



Z = j DMexv[-\ tr ( M^) - ^ tr ( M% 



The integration in ( p.85| ) is over A^xA^ hermit ian matrices. According the 't Hooft diagram 
technique the perturbative expansion in the coupling constant of the correlation functions 
( |2.85| ) is given by a sum of all planar double-line graphs Due to a normalization factor 
external lines corresponding to trM^" can be treated as lines of a generalized 
vertex. Note that all vertices of these graphs have the equal orientation (left or right 
in dependence on an adopted convention). We shall call two double-line planar graphs 
topologically equivalent if one of them can be transformed into the other by a continuous 
deformation on the plane. 

A planar non-vacuum graph is an half-planar graph if it can be drawn so that all 
its vertices lie on some plane line in the right of the generalized vertex trM^" and all 
propagators lie in the upper-half plane without overlapping pGf. Also we shall call a 



planar graph the half-planar irreducible graph if it can be represented as an half-planar 
graph in an unique way. A simple analysis shows that an half-planar graph is the half- 
planar-irreducible one if it does not contain any tadpole subgraphs. By a graph with a 
tadpole, as usual, we mean a graph with a subgraph which contains two lines coming 
from the same vertex so that after removing of these two lines the remaining subgraph 
becomes disconnected with the rest of the graph. 

In |2^ it has been shown that if one considers some special approximation for the 
planar theory, namely so called half-planar approximation, then Green's functions of the 
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one-matrix model in this approximation coincide with correlation functions in a corre- 
sponding Boltzmann theory, 

nf/(/=<0'"(l + #')-^>. (2.86) 

In both hand sides of (|2.86|) we omit all graphs with tadpoles. The raison to omit the 
graphs with tadpoles is that otherwise we have a double- counting for graphs with tadpoles 
for the Boltzmann theory. We call this approximation the half-planar approximation. 

A pure combinatorial proof of equation ( |2.86| ) is based on the following two statements. 
The first one states that the planar correlation functions without vacuum insertions in all 
order of perturbation theory are represented by a sum of all topologically non-equivalent 
graphs without any combinatoric factors. The proof of this statement is presented in 
Appendix B. According the second one in the Boltzmann theory all graphs contribute 
into correlations functions without any combinatoric factors. 

Let us compare the half-planar approximation with the planar approximation in first 
orders of perturbative expansion in the coupling constant. We collect corresponding 
diagrams on Figures 5-7. 

The explicit formula for an arbitrary planar Green's function without vacuum inser- 
tions is 134 



^2n{g) = .f"^' a2"(2n + 2 - na'), (2.87) 
n\[n + 2)1 

where 



a' = —[^l + l2g~l]. (2.88) 



For n = 1, 2, 3 we have 



n2 = ^a2(4 -a') = l-2g + 9g^ - 54/ + 378/ - ... , 

n4 = a^(3 -a'^) = 2-9g + 5Ag^ - 378/ + 2920/ - ... , 

Ue = a%8 - 3a^) = 5 - 36g + 270/ - 2164^^ ^ _ _ _ ^2.89) 

On Table 1 the results of numerical calculations of half-planar Green's functions Q2, 
^4, Qq and the planar Green's function 112, n4, Hg for some values of the coupling constant 
g are presented 



Table 1 



9 


10-3 




10-2 




10-1 


1 




10 




102 




103 




n2 


9.98- 


10-1 


9.81 ■ 


10-1 


8.6 ■ 10-1 


5.2- 


10-1 


2.1 ■ 


10-1 


7.4 ■ 


10-2 


2.4 ■ 


10-2 




9.97- 


10-1 


9.71 ■ 


10-1 


8.0-10-1 


4.0- 


10-1 


1.3 ■ 


10-1 


3.2 ■ 


10-2 


7.4 ■ 


10-3 




1.99 




1.96 




1.42 


4.84 


•10-1 


7.87 


■10-2 


9.26 


■ 10-3 


9.76 


■ 10-^ 




1.99 




1.91 




1.38 


4.43 


■10-1 


7.16 


■10-2 


8.65 


■ 10-3 


9.37 


■ 10-^ 


He 


4.96 




4.67 




2.92 


5.48 


■10-1 


3.47 


■10-2 


1.38 


■ 10-3 


4.70 


■ 10-5 


Ge 


4.97 




4.70 




3.11 


7.78 


■10-1 


9.64 


■10-2 


1.00 


■10-2 


1.00 


■ 10-3 



One can see that the answers for half-plane Green's functions G2, Ga practically satu- 
rate the planar Green's functions 112, n4 in board interval of the values of g. Ge is in the 
good accordance with Hg for g only for small g. 
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Let us make also a comparison of these two theories for the case when all tadpole 
graphs are neglected. The matrix theory the rejection of tadpole graphs is equivalent to 
adding to the Lagrangian the term 

-gIl2tiiM^), (2.90) 

where 112 is tadpole-free two-point planar Green's function. So the tadpole-free planar 
Green's functions are defined as 

n^. = Jhn ^ I DMtr{M'-) exp[-l(l - 2gfl2)tr{M') - ^tr{M% (2.91) 

The perturbative expansion in the theory (|2.91| ) is valid only for 

l-2^n2>0. (2.92) 
For this case making the change of variable 

M = (1 - 2gfl2)~^M 

in the (|2.91|) one gets 

Hsn = / DM{1 - 2gfl2)-''tr{M'nexp[~^tr{M') - -^tr{M% (2.93) 

where 

9 



(1 - 2gU2Y 

So the Green's functions Il2n{g) in that all tadpole contribution are dropped out can be 
expressed in terms of full Green's functions 'n.2n{g) 

Ii2n{9) = (1 - 2(;n2)-"n2„( I ). (2.94) 

(1 - 2^n2)2 

Now let us reformulate the Boltzmann theory in such way that it reproduces only 
tadpole-free half-planar graphs. One possible way to do this is to consider the formulation 
of Boltzmann theory using two fields. 

For this purpose let us defined the fields ip and as follows 

= a + a+, ^ = b + b^, (2.95) 
where a, a^, b and b^ satisfy the following relations 

aa^ = 1, bb^ = 1, 

ab+ = 0, ba+ = 0. (2.96) 

This algebra has the realization in the Boltzmannian Fock space. For two degrees of 
freedom the vacuum |0) satisfies 

a|0) = 0, &|0)=0. (2.97) 
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Now let us define an interaction Green's functions by tlie formula 



J- "T Oint 



|0), 



(2.98) 



where the interaction Sint has the form 



Sint = gi) : (jxp : i) = g^cpcp^p - gi>i>. 



(2.99) 



One can see that the connected tadpole-free parts of Green's functions in Boltzmann 
theory described in subsection 2.2 coincide with the connected parts of Green's functions 

Id 
,1 

F„ = = (01^ : : ^ . ^ 10). (2.100) 



1^2— field 
t^l,n-2,l 



J- + ^int 

Let us write down the system of Schwinger- Dyson equations for the Green's functions 
F2 and F4 

F2 = 1 - gF^, 

= - gFl (2.101) 

The results of numerical calculation of F2 and F4 related to a perturbative branch of 
solution of ( |2.101|) are presented on Table 2. This solution has a phase transition point 
5^0 = 0. When g goes to go then F2 and F4 go to infinity. Also on Table 2 we present 
the numerical values for 62 and nconnected ^ jj^ _ fQj, g ^ = 0.421875, where go is 
determined from the following equation 



1 - 2^n2 = 0. 



(2.102) 



Table 2 



9 


10-3 




10-2 




10-1 


0.4 


0.421870 


0.9 


0.99 




1.00 




1.00 




1.01 


1.14 


1.16 


2.82 


9.56 


n2 


1.00 




1.00 




1.01 


1.16 


1.18 






Fj 


-9.99 


10-4 


-9.90 


10-3 


-0.0925 


-0.359 


-0.382 


-2.03 


-8.64 


ytconnected 
114 


-9.98 


10-4 


-9.81 


10-3 


-0.0872 


-0.394 


-0.439 







Perturbation series for F2, F4, II2 and n™""'^'^*'^'^ have the forms 



-l+g + ^l + 2g-3g'' 



l + g-2g^-Vl + 2g-3g'' 
2g\l-g) 



l + g^-g^ + 3g'~6g'' + ... 



-g + 2/ - 5^3 + _ ^5^5 ^ 



(2.103) 



Ha = 1 + / - 2^3 ^ - 42/ + ... 



connected 



-g + 2g^ - 10^3 ^ 42^4 _ 299^5 _^ 



(2.104) 
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3 Boltzmann Theory for Finite Number of Degree of 
Freedom 



3.1 The Schwinger-Dyson Equations 

In this section we will study the Boltzmannian Schwinger-Dyson equations for the theory 
with finite number of degrees of freedom. 

First of all let us derive the Boltzmannian Schwinger-Dyson equations. This can 
be done similarly to the zero-dimensional case. We adopt the following notations. Let 
a,b,c, ... are multi-indices which label degrees of freedom, 0(a) = (p'^ia) + (f)~{a), where 
4>^{a) and (f)~{a) satisfy the relation 

0-(a)0+(6) = G(°)(a,6) (3.1) 

and G^^\a, h) is a matrix which has the inverse K{a^ b): 

J2K{a,c)G{c,b) =6{a,b). (3.2) 

c 

The algebra (|3.1|) has the realization in the Bolzmannian Fock space generated by the 
vacuum |0), 0~(a)|O) = 0, and n-particle states 

|ai,...,a„) = 0+(ai)...0+(a„)|O). (3.3) 

Note that states ( |3.3| ) have no any symmetry properties under permutation of particles 
because there is no relations between 0"^(ci) (f)^{b) and (f)'^{b) 0+(a). 
Free correlation functions 

G(°)(«i,-,«n) = (O|0(ai)...0(aJ|O) 

in the Bolzmannian Fock space satisfy the relations 

c 

n 

= J2 '5(afc,ap)G(°)(afc+i,-,«p-i)G'(°n«i,- n ) + 

p=k+l 
fc-1 

6{ap,ak)G^°\ap+i, ...,ak-i)G^°\ai, ap_i, a^+i, ...,a„), (3.4) 

p=i 

that follow from the Wick theorem in the Boltzmannian Fock space. We also call equations 
( |3.4]) the free Schwinger-Dyson equations. 

Interacting Green's functions Q{ai, an) are defined by an analogy with the 0-dimen- 
sional model: 

C;(ai,...,a„) = (O|0(ai)...0(a„)^^(0)|O), ^](0) = ] (3.5) 

Let us note that the Green's functions ( p.5| ) have a symmetry under the permutation of 
arguments in inverse order 

Q{ai, an) = G{an, ai)- (3.6) 
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There are no another symmetry properties. 

The form of Vint{(l>) depends on the model under consideration. We are restricted 
ourself to the case with Vint in the form 

Vinti^)=gY^i<Pia)r. (3.7) 



To derive the Sc winger- Dyson equations for correlation functions ( |3.5| ) let us rewrite 
the Green's function ( p.5|) in terms of the free Green's functions, 

oo 

g{a^,...an)=Yl E (-^?r(O|0(ai)...0(a.)(0(&i)r-(0(Ur|O) = 

m=0 bi,...bm 

oo 

G^^\ai, ...,an,bi, ...,bi, ...,bm, ...,bm) (3.^ 

m=Obi,...bm 



Using the free Schwinger-Dyson equations (|3.4| ) we have 



Y K{ak, c)Q{ai, a^-i, c, Ofc+i, a^) = 

c 

oo 

Y K{ak, c)Y E (-fi')™'^'"°''(«b flfc-i, c, Ofc+i, an, bi, 61, 6™, bm) = 
c m=obi,...,bm ' ' 

= Ji + Jn + Jm, (3.9) 

where 

00 fc— 1 

Ji=Y E E(~fi')"'^(S''^fc)G'(°)(ai, ...,ap_i,afc+i, ...,a„,6i, ...,61, ...,6„)x 
m=ofei,...,6„p=i ' ^: ' ' ' 

xG'(°)(ap+i,...,afc_i) (3.10) 



Jn=Y E E (-5')™'5(afc,ap)G'-°n«i' •••''^"'^1' •••'^1' •••'^m' •••'^m)^ 

m=o p=fc+i j:; ' ' 

xG'(°)(afc+i,...,ap_i) (3.11) 

m r 

Jill = Y J2Yj(~9)'^^(^k,bi)G^^\ai, ...,afc-i,6i, ...,61, 6;, ...,bi)x 

bi,...,bm 1=1 s=l " ' ' ' 

X G^°\bi, ...,bi,ak+i, ...,a„,6z+i, ...,6^, •••,&m)- (3-12) 

Using the relation 



m=0 1=1 1=1 m=l 
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and making the summation over 6/ and then changing the summation index 

m ^ m' = m — I, 
one can represent the term Jm in the form 

oo oo r 
1=1 m'=0 bi s=l 



s-1 



- gY] S{ai, Ofc-i, afc, Ofc, )^(afc+i, a„) (3.13) 

s— 1 

Putting all together we obtain the half-plane Schwinger- Dyson equations 

c 

fc-1 

X 5(ap, afc)^(ai, ap_i, a^+i, a„)G'^°^(ap+i, afc_i) + 

n 



p=fc+i 



fl'H^loi) •••,afc)^(afc, ...,0^,0^+1, ...,a„). (3.14) 



s=l 



Also we have the equations 

Q{ai, ...,an) = G^^\ai, ...,a„) - gJ2^i^^^ ■■■,an,b, ...,b) (3.15) 

which follow from the average the operator identity 

1 - Vint 



1 + Vint 1 + Vint 



1 (3.16) 



Let us consider the interaction g J2a[4>{(^)V^ write down the Scwinger-Dyson equa- 
tions ( |3.14|) corresponding to (n, k) = (2, 1), (4, 2). We have 

Y K{a, e)Q{e, b) = 5{a, b)Q — gQ{a, b)Q{a, a) — gQ{a, a, a, b)G, (3-17) 

e 

Y K{b, e)Q{a, e, c, d) = 6{b, c)Q{a, d) + 5(a, 6)^(c, d) — 

e 

gQ{c, d)g{a, b, 6, b) - gg{a, b)g{b, 6, c, rf), (3.18) 

where Q is vacuum Green's function. Equation ( |3.17| ) and ( |3.18| ) together with the equa- 
tion for vacuum Green's function that follows from (|3.15|) 

^ = 1 — 5? X ^(a, a, a, a) (3.19) 

a 

form the closed system for Q, Q{a, b) and Q{a, b, c, d). 
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3.2 Diagonal Kinetic Operator 

In this subsection we will study the theory with K degrees of freedom and diagonal kinetic 
operator. We assume the following form of Lagrangian 

^=\j:^i+9j:<- (3.20) 

^ a=l a=l 

The latin indices a, 6, c, ... run the values 1,2, ...,K. The operators are the sum of 
Boltzmann creation and annihilation operators = $^ + which satisfy the relation 

K<^t = Sa,b. (3.21) 

The Green's functions are defined by 

g(a,&,...,c)=g.b...e=(0|$a$fc-$c ^ J O)- (3.22) 

^ ' 9 2-ia ^ a 

It is evident that 

Ga^ = Gi)j) (3.23) 

" n 

for any a and b and 

Qabb...b = (3.24) 

for a b. 

Let us write down the equations (|3.17|) , (|3.18|) and ( p.l9|) for this case: 

Gaa = G - gQla " QQQaaaa, (3.25) 

(3.26) 

g=l- KgGaaaa- (3.27) 

From these equations one can write the equations for Q: 

4gg^ + 12{K-l)gg'^ + [12{K-lfg + K]g^ + [2K{K-l)+Ag{K-lf]g-K{2K-l)=0. 

(3.28) 

This equation can be solved and the Green's functions and gaaaa are expressed in 
terms of g. 

3.3 Two Degree of Freedom and Non-Diagonal Quadratic Ki- 
netic Operator 

In this subsection we shall consider the theory with two degree of freedom and non- 
diagonal kinetic operator. The Lagrangian has the form 

C = WK,,t^a<^t + 9Y.^t- (3.29) 

a,b a 

The indices a,b,c, ... will take the values 1 and 2. We adopt the following parametrization 
of the kinetic operator Kab 

Ku = K22 = a, Ki2 = K2i = -p. (3.30) 
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In this case we have 



$1 $+ = $2 $2 



a2 -/?2 



$1 ^\ = $2 $1 



(3.31) 



The Lagrangian is invariant under the transformations $i $2) ^2 ^ ^i- Therefore 
all Green's functions are invariant under the change of indices 1 — 2, 2 — 1. In particular, 
we have 

Gab = Gba, (3.32) 

Ga...a = Gjj J) (3.33) 



for any a, b and n. 

Let us write (|3l9D , §J7\) for (a, 6) = (1, 1), (1, 2) and ( glSD for (a, 6, c, d) = (1, 1, 1, 1), 
(2,1,1,1), (1,2,1,1), (2,2,1,1): 



G = 1 — gGnn — gG2222, 



(3.34) 



aGii - PG21 = G - gGii - gGiuiG, 
aGi2 - (^G22 = -gGuGu - gGuuG, 

aGnn — /?^i2ii = '^Gn — ^gGuGiin, 

«^2iii - /5^22ii = G21 - gGiiG2iii - gG2iGiiii, 

— PGllU + «^1211 = — 5'^11^1222 — 5'^12^2211, 
— /5^2111 + ttfell = Gil — gGllG2222 — fl'^22^2211- 

The following symmetry properties of Green's functions are evident: 

Gu = G221 G12 = G211 Gun = ^2222, 



(3.35) 
(3.36) 
(3.37) 
(3.38) 
(3.39) 
(3.40) 



^1211 — ^2122; ^1122 — ^2211; ^1112 — ^2221 — ^1222 — ^2111- (3-41) 

So the system ( p.34|) -( p.40|) contains 7 independent equations and 7 unknowns variables 
G, Gu, Gi2, Gnu, ^1112, ^1211, ^1122 and therefore can be solved. 

To get the equations for Green's functions without vacuum insertions one has to set 
^ = 1 in (|3.35| ) - (13.401) . Then we have a system of 6 equations for 6 variables. The 
numerical solutions of this system for a = 1 and some values of (3 and g are presented on 
Tables 3-5. 



Table 3. (3 = 0.01 
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9 


iU 










iU 


i 




1 n 




1 n2 


1 n3 
iU 








Q Q7 




Q 70 

y. / z • 


1 n~ 


1 


7 QQ 1 n— 1 
/ .yy • lu 


o.yo • 




l.ZD • 




Q 9 1 n-2 

O.Z • lU 


7 


1 H" 
lU 


-3 




9.94 


•10-3 


9.43 • 


10- 


3 


6.38 • 10-3 


1.58 • 


10-3 


1.58 • 


10-^ 


1.03 • 10-5 


1.51 


10- 


-6 


Ciiii 


1.99 




1.91 






1.38 


4.43 • 


10-1 


7.16 • 


10-2 


8.65 • 10-3 


9.37 


10" 


-4 


G^1112 


1.98 
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The dependence of the 2-point Green's functions Gn and G12 on the coupling constant 
g for different values of /3 is also presented on Fig. 8. 

The limit /? ^ corresponds to decoupling of the degrees of freedom. In this limit 
the mixed Green's function G12 goes to zero and we reproduce the case of one degree of 
freedom. For (3 = 0.9 the mixed Green's function is of the same order as the diagonal 
Green's function. 

4 Boltzmann Correlation Functions in D-Dimensio- 
nal Space-Time 

4.1 The Schwinger- Dyson Equations 

In this section we will derive the Schwinger-Dyson equations for Boltzmann correlation 
functions in L>-dimensional Euclidean space. 
We adopt the following notations. Let 

(j){x) = (p'^{x) + (p-{x) (4.1) 
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be the Bolzmann field with creation and annihilation operators satisfying the relation 



0-(x)0+(i/) = D(x,i/), (4.2) 

where 

D{x, y) = — / rf^P^— ^ (4.3) 

is D-dimensional Euclidean propagator. In what follows we will assume that 0(a;) is 
a scalar field. A generalization to the case when carries isotopic or tensor indices is 
straightforward. 

By analogy with the case of finite number degrees of freedom we consider the following 
correlation functions 

Fn{Xn,...,Xi) = (O|0(Xn)...0(Xi) _^ . \ , .J O), (4.4) 

where Cint{<P{x)) is a local invariant polynomial in the field and its derivatives. As before 
we restrict ourself by the case of the quartic interaction £j„t(0) = gcf)^. 

Using the results of the previous section one can immediately write down the Boltz- 
mannian Schwinger- Dyson equations for correlation functions ([4.4|) . We have 



[- A +m^)^2,F2„(x2n, ■■■Xl) = -g{F2l{x2U X2I-I, ■■■Xl)F2n-2l+2{x2n, ■ ■ ■X2I+I, X2I, X2l) + 

F2l+2{X21, X2hX2U X2I-I, ...Xi)F2n-2l{x2n, ■■■X2l+l)) + 
S{X21 - Xi)F2n-2l+i-l{x2n, ■ ■ .X2I+I, Xi-i, . . .Xi) F^fli_^{x2l-l, ...Xi+i) + 



i<2l 



6{X21 - Xi)F2n+2l-i-l{x2n, ■■■Xi+1, X2I-I, . . .Xi )F>_^2/-1 (^^i-l 5 ■■■X2I+1), I < n, (4.5) 

2l<i 

and the similar equations for X21+1 

(- A +m^)x^^^^F2n{x2n, ■■■Xl) = 

-g{F2l{x2l, X2I-I, ...Xi)F2n-2l+2{x2n, ■■■X21+2, X2I+I, X2I+I, X2l+l) + 

-^2«+2(a^2«+l5 X2I+I, X2I, ■■■Xi)F2n-2l{x2n, ■■■X2I+2, a;2/+l)) + 

Y S{X21 - Xi)F2n-2l+i-2{x2n, ■■■X21+2, Xi-l, . . .Xi) F^ili{x2l , ...Xi+i) + 
i<2l+l 

Y ^{^21 - Xi)F2n+2l-i{x2n, ...Xi+i,X2l, . . .Xi)F/°),,„2 (^^i-l ) ■■■X2I+2), I < U. (4.6) 
i>2l+l 

Let us compare this set of equations with the usual Schwinger-Dyson equations and 
with the planar Schwinger-Dyson equations. The usual Schwinger-Dyson equations for 
scalar field with ip^ interaction have the form 

(- A +m^)^^Gn{Xn, --Xl) = -gGn+2iXn, Xn, Xn, Xn-1, ...,Xi) + 
n-1 

Y ^{Xn - Xi)Gn-2{Xn-l, . . .Xi-i, Xi+i, Xl) (4.7) 
i=l 



25 



and Gn{xn, ...xi) is a symmetric function of its arguments. 

To write down the planar Schwinger-Dyson equations one has to modify in ( |4.7|) only 
the Schwinger's terms: 

( ^ -\-ni )x'„nn('^n; Xfi, Xij gYln^2(^-^m -^nj -^ny -^l)"!" 

n-1 

+ ^5{xn- Xi)ni_i(xi_i, xi)n„_i_i(x„_i, ajj+i), (4.8) 

where 

n„(a;„, ...xi) = Jiin < tr ( M{xn)...M{xr)) > (4.9) 

and Iln{xn, ...Xi) is invariant under cyclic permutations of its arguments. The Schwinger 
terms for the Boltzmannian Schwinger-Dyson equations differ from the Schwinger terms 
for the planar Schwinger-Dyson equations. From ( |4.5| ) we see that the Schwinger terms 
are linear with respect to the interacting Boltzmann correlation functions meanwhile 
the Schwinger terms in the planar equations are quadratic with respect to the planar 
correlation functions. There is also a modification in the term representing the interaction. 
Instead of the linear term for full or planar equations we have quadratic terms in the 
Boltzmannian Schwinger-Dyson equations. The usual Schwinger-Dyson equations may 
be written in the operator form as follows 

5{So + S,rrtM ^ Schwinger's terms = 0. (4.10) 

0(f{x) 

The planar Schwinger-Dyson equations in the operator form can be presented in terms 



of the master field $ and the master momentum 11 0, |T9|, |T8|, ^ 



S^{x) 

where S is an action and $ and 11 satisfy [El p7 



+ 2n(a;)]|0) = 0, (4.11) 



[U{x),^y)]=t6^''\x-y)\0){0\. (4.12) 
The Boltzmannian Schwinger-Dyson equations can be written down schematically as 

1>{x)\n){0\ f'"'^^^ + l»>(0| '^f!"'''^' + Schwinger's terms = 0. (4.13) 

We see that the Boltzmannian Schwinger-Dyson equations in the operator form are more 
complicated than the usual Schwinger-Dyson equations since the interaction terms contain 
projectors on the vacuum state. But as we will see in the next subsection just due to 
these projectors it will be possible to get a closed set of integral equations for two- and 
four-point Green's functions. 
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4.2 A Closed System of Equations for Two- and Four-Point 
Green's Functions 



Let us write down equations ( |4.5| ) and ( |4.6| ) for the two-point correlation function 

(- A +m'^)^F2ix,y) = -gF2{x,y)F2{x, x) - gFi{x,x,x,y) + 5{x - y), (4.14) 



A +m^)yF2{x, y) = -gF2{y, y)F2{x, y) - gFi{x, y, y,y) + 5{x - y), 



(4.15) 



and the four-point correlation function 

(- A +m^)^Fi{x, y, z, t) = -gFi{x, y, z, t)F2{x, x) - gFQ{x, x, x, y, z, t) + 

S{x - y)F2{z, t) + 5{x - t)F!f\y, z), (4.16) 

(- A \rr?)yF^{x, y, z, t) = -gF2{z, t)Fi{x, y, y, y) - gF^iy, y, z, t)F2{x, y) + 

6{x - y)F2{z, t) + 6{y - z)F2{x, t), (4.17) 



(- A +m'^)^F^{x,y,z,t) = -gF2{z,t)Fi{x,y, z, z) - gF^{z, z, z,t)F2{x,y) + 
6{z-t)F2{x,y) + 6{y-z)F2{x,t), 



(4.18) 



(- A +m'^)tFi{x, y, z, t) = -gF^ix, y, z, t, t, t) - gF2{t, t)Fi{x, y, z, t)) + 



5{z-t)F2{x,y) + 5{x-t)F!f\y,z). 



(4.19) 



Here we assume that all vacuum insertions are dropped out. Note that not all these 
equations are independent. Indeed, having in mind the symmetry of half-plane Green's 
functions under the permutation of their arguments in the inverse order ( p. 61 ) one can 
reproduce equations ([4.15|) , ( [4.17|) and (|4.19| ) from equations ( |4.14|) , ( [4.18|) and [4.16| ), 
respectively. 

We see that equations (|4.18|) and ([4.17|) in accordance with the general formula ( [4.13|) 
do not contain six-point correlation functions. Let us show that from ( ^4.17 ) and ( [4.14 ) 
one can get a closed set of integral equations which is sufficient to find F2 and F^. To this 
purpose let us rewrite ( [4.17|) in the integral form 



F,(., y, = -9 / fl(y, .)|F.(., t)F.(-, + f.K 

D(x, y)F.,(z, t) + D(y, z)F.,(x, t) (4.20) 

(see Fig. 10). The integral representation for the Schwinger- Dyson equation (|4.14|) is 
drawn on Fig. 9. As has been mentioned above the Boltzmann correlation functions 
do not possess any symmetry property except the invariance ( p. 61) . Therefore drawing 
correlations functions one has to mark the order of the external lines. The possibility to 
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do this is to draw arguments of the correlation function on a hne as it is done on Figs. 9 
and 10. 

For the two-point correlation function we have 

F2{x, y) = -g J[D{x, u)F2{u, u)F2{u, y) + D{x, u)Fi{u, u, u, y)]du + D{x, y), (4.21) 
(see Fig. 9). 

From (|4.2CI|) and (|4.21|) one gets the relation between J-'i{x, y, z, t), the connected part 
of -^4(0;, y, z, t), and F2{x, y). To write down this relation note that 

F4(x, y, z, t) = T^[x, y, z, t) + F2{x, y)F2{z, t) + F2{x, t)D{y, z), (4.22) 

(see Fig. 11). Finally we get 

J^4{x,y,z,t) = -9 J F2{x,u)D{y,u)[D{u, z)F2{u,t) + J^^{u,u, z,t)]du (4.23) 

F2{x,y) = D{x,y) - g J D{x,u)[D{u,u)F2{u,y) + 2F2{u,u)F2{u,y) + J^i{u,u,u,y)]du. 

(4.24) 

Let us write down this system in a more suitable form. Note that there are tadpole 
contributions in the right hand side of equation ( [4.24| ). It is convenient to include these 



terms into a mass renormalization. To this purpose let us substitute the representation 
(lO^ ) in ( ^:T5D . We have 



(- A +m^)yF2{x,y) = -g{2ci + co)F2{x,y) - gj^i{x,y,y,y) + 5{x - y), (4.25) 
where ci = F2{x,x), cq = D{x,x). Including these constants in a new mass M we get 

F2{x, y) = Dm{x -y) -g J duJ^^{x, u, u, u)Dm{u, y). (4.26) 

Instead of equation ( [4.23| ) it is more convenient to deal with an equation for a partially 
one-particle irreducible (IPI) 4-point function r4(x, t) which is defined as 

14 = F^'D-'D-'F^^T,, (4.27) 

or more precisely 

T^{x,y,z,t) = J F2\x,x')D-\y,y')D-\z, z'))F^\t,t')J^4{x',y',z',t')dx'dy'dz'dt'. 

(4.28) 

Note that in the right hand side of ( [4.28 ) we multiply JF4 only by two full 2-point Green's 



functions while in the usual case to get IPI Green's function one multiplies an n-point 
Green's function by n full 2-point functions. From ([4.23|) and ( [4.28|) we get 



Ti{p, k,r) = -g-g J dk'F2{p + k- k')D{k')Ti{p + k - k' , k' , r). (4.29) 
Equation ( [4.26|) can be also rewritten in terms of r4 



F2{p) = Dp-gJ dk'dk"D{k')D{k")F2{p - k' - k")r4{p - k' - k" , k' , k"). (4.30) 
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Equations ( |4.29| ) and ( [4. 301 ) are presented graphically on Fig. 12. 

Equation ( |4.29|) is the Bethe-Salpeter-like equation with the kernel which contains an 
unknown function F2. As in the usual case we can write down F2 in term of self-energy 

and write down equation ( [4. 261 ) as an equation for S2, i.e. 



^2{p) = J dkdqF2{k)D{q)D{p -k- q)T,{p, k, q). 



(4.32) 



Equation ([4.32|) is similar to a usual relation between the self-energy function S2 and 
4-point vertex function for field theory. Note that a closed form of the representation 
for 4-point vertex function ( 4.29 ) is a distinctive feature of the Boltzmann field theory. 

One can choose an approximation in the system of equations ( |4.32D , ( f4.29| ) which gives 
a so-called rainbow approximation. To get this approximation one takes r4 ^ —g and 
restricts by the first two terms in the relation ([4.31|) , 



pappr ' 



M2 



(4.33) 



Within this approximation the solution Sa^^^ = 11 of equations ([4.32|) is reduced to the 
solution of the following equation for 11 



n{q) = Y[,{q)+g^ j 



d\' I{q - q') 



rn(g') 



(27r)4 (g2 + M2)2 

where no(g) denotes the second order contribution to the self-energy 



(4.34) 



d'^ki f d'^k2 



and /(g) is a "potential" 

J(p) 



1 



(27r)4 {kl + w?) {kl + m2)((g -ki + kif + M2) ' 



(4.35) 



d\ 



(27r)4 (g2 + m?){{p - g)2 + m?) 



(4.36) 



Equation (|4.34| ) is known as the rainbow equation and it was studied by Rothe p2 



The integrals ( ^.36| ) and ( |4.35| ) contain divergences. Performing a subtraction one gets 
the following expression for the renormalized "potential" (for simplicity we assume that 
m = 0) 



1 1 



(4vr) 



(4.37) 



with K being an arbitrary subtraction point. The renormalized second order self-energy 
has the form (eq. (A. 9) from ||3^ ) 



no(g) 



2(47r) 



{A + Be + (f - ^) Ml - ^) + 2| ^1 1„(1 - x-^)}. (4.38 



where A and B are arbitrary renormalization constants. 
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For the case of exact equations ( [4.29| ), ( |4.31| ) and ( |4.32| ) one also has to perform 
appropriate renormalizations. There are two possibilities to discuss renormalizations. 
One can do an appropriate subtractions directly in the integral equations or one can do 
them order by order in the perturbative expansion. Both ways are equivalent but the 
second one is more closed to renormalizations of usual field theories. We are going to 
discuss it in the next subsection. 



4.3 Renormalization of the Boltzmann Field Theory 

In this subsection we will describe the renormalization of the Boltzmann field theory. To 
make a formal expression for correlation functions 

f„(x.....x„) = (0|^(^i)-^(^-.) ^^^,„^^^^^^(^(^)) |0>, (4.39) 



finite we apply R-operation |33]. In ( 4.39| ) has the decomposition ( [4.1|) with (j)~^{x) 



and (f)~{x) satisfying ( |4.2| ) and the interaction Lagrangian is an invariant polynomial of 0"'', 
(f)~ and its derivatives. The perturbative expansion of Boltzmann correlation functions 
(|4.61| ) is represented by the sum of half-planar Feynman graphs {H}. We write Fh for 
the unrenormalized value of graph H 

Fn^Pu-.-Pn) = j dki...dkilH{ki, ...,ki;pi, ...,pn), (4.40) 

where pi, ...,Pn are external momenta and the integration is over internal momenta ki, ki. 

Let us apply the standard i?-operation to Boltzmannian Green's functions (|4.39| ). If 
H is not a renormalization part {H has no overall divergence) we set 

Rh = Rh, (4.41) 

and if is a renormalization part {H is IPI and has an overall divergence) 

i?H = (1 - Th)Rh, (4.42) 

where Th is a subtracting operation which fixes the renormalization scheme. The function 
Rh is defined recursively by 

Rh = Ih+ E lH/{,„...,,.}fl{-TM. (4.43) 



The sum in ( 4.43 ) is extends over all sets {71, 7s} of renormalization parts of H which 
are mutually disjoint and different from H. The reduced diagram H/ {71, 7s} is defined 
by contracting each 7j to a point. The explicit formula for the renormalized integrand is 
given by the forest formula, 

RH = il + j:Ui-T,))lH (4.44) 

u -yeu 

with the sum going over all iY-forests. Up to this point a specific character of the Boltz- 
mann field theory has not been taken into account. 
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To prove the renormalizability for some given Cint we have to prove that all sub- 
tractions can be collected in corresponding counterterms. At this point we meet with a 
specific character of the Boltzmann field theory. There is the following subtlety in this 
step. Half-planar graphs contain renormalization parts which are half-planar as well as 
renormalization parts which are not belong to the set of half-planar graphs. With an 
half-planar renormalization part we have to associate the counterterm 

£„(i/„) = {-TnAJr- (4.45) 

Note that there is no the combinatoric factor. We will see that with non-half-planar 
renormalization parts we have to associate counterterms which are not functions of but 
are functions of 0^. 

To specify our discussion let us consider correlation functions ( [4.61| ) for the cf)^ inter- 



action in D = 4. To avoid problems with tadpoles we consider the two-field formulation 
where 

ip{x) = ip~ {x) + ip^i^x), 0(x) = 0~(a;) + (p'^^x), 
ijj- {x)ip+ (y) = DMix,y), 0"(x)0+(?/) = D^{x,y), 
V^^(x)0+(y) = 0"(x)^+(?/) = 0, 

DM{x,y) = J ^-^.^'''^-''^ ^ = iM,m). (4.47) 

We draw 0-propagator by thin lines and '0-propagator by thick lines. Only two- and 
four-point graphs may be divergent. Any two-point subgraph of the half-planar graph 
is also from the set of half-planar graphs and performing contractions of the two-point 
subgraphs of the half-planar graph one gets again a half-planar graph. But four-point 
subgraphs of a half-planar graph may be of half-planar type and may be not. We refer 
to the later CclSG clS to the case of subgraphs of Il-type. Particular examples of divergent 
parts of Il-type of the two-point half-planar graph as well as an examples of half-planar 
divergent subgraphs are presented on Fig. 13. 

At the first sight a four-point divergent part of Il-type requires a counterterms ■0^- 
But a careful examination of half-planar diagram shows that in fact only the counterterms 

: ij^j : (4.48) 

are necessary. Indeed, let us consider a graph with a four point divergent part of Il-type. 
There are two possibilities to make a contraction of this graph to a point, namely, making 
a contraction to the left vertex A, or making a contraction to the right vertex B (see 
Fig. 14). In the both cases after the contractions we are left with half-planar graphs with 
the insertions 

: ^0 : ^'iIj- (4.49) 

and 

^+0+ : :, (4.50) 
respectively. One has to choose one of these possibilities. Let us take the first one. 
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The above consideration shows that reduced graphs if/ {71, ...7s} are always from the 
set of half-planar graphs. This fact permits to collect all divergent parts to counterterms 

£,^(^+,^-,0) = (4.51) 

H 

where the summation is over all renormalization parts. Therefore 

= {Z^ - l){dij^f + 6Mi;^ + 6 : <j)<f) : i; + 6 X : ^j^j : ^j-^-. (4.52) 



We add counterterms ( [4.52| ) to the interaction Lagrangian and consider the renormal- 
ized correlation functions 

(O|^(Xi)0(x2)...0(Xn-i)V^(Xn) , . ,^ r , \ , M + / - ' ^^■^^'> 

1 + J d^xlgi/j : (j)(f) : ijj + Cct{w, ^ , 0)J 

It is easy to see that the counterterms remove all divergencies in the graphs corre- 
sponding to (|4.53|). Indeed, let us consider some half-planar graph If n > 4 then Hn 



is IP reducible and it does not require any subtraction. If H2 and Hj^ are overall divergent 
graphs then these divergences are subtracted by the graph created by the counterterm 
vertices {Z^ — l){dip^y, SMip"^ and dgip : (fxp : ip respectively. Now let us supposed that 
Hn has overall divergent IPI subgraphs 71,. ..,7^. Since the graphs 71, ...,7^ are overall 
divergent there are corresponding counterterms in ( [4.51j ). Divergent subgraphs may be 



halp-planar ones, that corresponds to the first three terms in ( [4.52|) , or Il-type graphs. 



that corresponds to the last counterterm (|4.52|). These counterterms produce the graphs 



in which one or more of overall divergent subgraphs are replaced by their counterterms 
and there are the corresponding terms in the forest formula. The counterterms ([4. 521) do 



not generated new divergences which are not cancelled with divergences of graphs gener- 
ated by the basic Lagrangian since there is one to one correspondence between the terms 
of the forest formula and the terms generated by the counterterms ( [4.52[ ). 



The distinctive feature of the Boltzmann theory is that there are the wave function 
counterterms {6Z — l){d'ilj)'^ and the mass counterterms Sm^tp"^ only for ip-field (see Fig. 
15). 

To associate counterterm ( ^.49| ) with a change of parameter of the theory it is conve- 
nient to add the corresponding term to the interaction Lagrangian, 

Cint{ip+, ip- , 4>) = gip : (p^ : ip + X : ipiJ : ij-ip'. (4.54) 



Therefore the theory (|4.54| ) is specifyed by four parameters, two masses M and m and 



two coupling constants g and A. 

Let us for completeness write down the Boltzmann Schwinger-Dyson equations for the 
interaction (4.54) 



(- A +M'^)^F2{x, y) = gF2{x, y)DM{x, x) - gF^{x, x, x, y)- 

XF2ix, y)F2{x, x) + AF2(x, y)DM{x, x) + 5{x - y), (4.55) 
(- A +m^)yF4{x, y, z, t) = -gF2{x, y)Fi{y, y, z, t) + 6{y - z)F2{x, t). (4.56) 
We see that equations (^4.56|) is just equation ( [4.23|) for the connected four-point Green's 



function. However equation ( [4.55| ) for the two-point function contains new terms. From 
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these equations it follows that equations ( |4.29| ) and ( ^.321 ) are not changed and instead of 
we have 

where 

Ml = M^ + A(F2(0) - Da/(0)). (4.58) 

Let us assume dimensional regularization and the minimal subtraction scheme with a 
scale fi. We have 

Cint = : 00 : 7/; + /i"-^A : : V^"^", (4.59) 

Cct = {Z^ - l){d^i^? + SM^i^^ + ^i^-'^dgij : 00 : 7/; + /i^-^(5A : i^^J : t/-"^", (4.60) 

F„(xi,...,x„;m, M,5(,A,/i) = (O|^(xi)0(x2)...0(x„_i)?/'(x„)(l + ^ (i^x(£i„t + /:rf))""^|0). 

(4.61) 

Now we can include the wave function renormalization Z^ — l'm the rescaling of the field 
%l) and the renormalized correlations functions up to the renormalization factor are equal 
to unrenormalized correlation functions in the theory with bare parameters, 

F„(xi, x„; m, M, g, A, /i) = Z^iF™'''^"(xi, m, Mq, ^o, Aq) = 

1 

+ / d°xCint{m, Mo, go, Aq 

where 



Z^^(O|^(xi)0(x2)...0(x._i)V^(x„) ; — 10), (4.62) 

^ 1 + d^xCir,Am, Mn, On, An 



Mo = Z^^(M2 + 5M2), ^7o = ^^V'"''(^7 + ^^7), Ao = Z//-^(A + a). (4.63) 

It is evident that the renormalized correlations functions being multiplied on do 
not depend on fi since F^"'^^"- depends only on bare parameters Mq,, go and Aq which 
do not depend on the scale parameter specifying the renormalization prescription. This 
implies the following renormalization group equation 

(/^^ + (^9-^ + /^A^ - ^m-q]^Jp + l)Fn{xi, ...,Xn] vfi, M, g. A, \i) = 0, (4.64) 

where 

The difference between ( [4.64|) and the usual renormalization group equation is that 



the anomalous dimension 7 in (4.64) is not multiplied on n/2. There is also a difference 



in the expression of the P and 7 functions in terms of the counterterms Sc and Sg. In 
dimensional regularization the counterterms are poles in [D — 4) 

As in the usual case one can get the low-order expression for the /^-function in terms 
of Ci and ai. The difference with the usual case comes from the fact that now there is a 
new expression for go in terms of 6g and Z^. We have 

go = fi'^''ig + 6g)Z;f;\ Aq = ^^'^(A + 5A)Z^2. (4.67) 
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Substituting ( |4.66| ) in ( [4. 671 ) gives 



go = [g -\ — — + poles ot order more them /J, 

X^ = ^^-D[x + M3i^ltz2^i9^ + polesof order more them 2]. (4.68) 

' dfj. 



Applying fJ^j- to both sides of ([4.68|) yields 



d 

Pgig) = (1 - Q\^^^9Ci{g, A) - ai{g, A)), 

/3a(A) = (1 - ^)(2Aci(^7,A) -6i((7,A)). (4.69) 
Taking into account the explicit form of Ci{g, A), ai{g, A) and bi{g, A) we get 

f^^^3,>^) = £-,+T^ + 0{g\X'), P.{gA) = ^,+0{g\\% (4.70) 



We see that the sign of the half-planar beta function (3g is the same as the sign of the 
standard beta function in theory but there is a numerical difference of the half-planar 
beta function with the standard beta function. 

A non-standard wave function renormalization in the Boltzmann theory ([4.4UD implies 



that a dependence of the renormalized correlation function on the scale parameter /x is 
given by the formula 

Fn{xi, ...,Xn; m, M{fi),g{fj.), A(/i), /i) = exp[- / —-f{g{fi"), A(/i"))] x 

J fi' jj, 

F„(xi, ...,Xn; m, M{fi"),g{fi"), A(/i"), /i")- (4.71) 

Since due to a special form of A-interaction A-vertices contribute only in the mass 
renormalization of ■i/'-lines and it is natural to expect that the dependence on A in the high- 
energy asymptotic of correlation functions may be neglected together with dependence 
on M(/i) and one get the following asymptotic formula 

F„(/€Xi, KXn, m, M{fi),g{fi), A(/i), /i) ~ /t''^~^(/t/i)F„(xi, Xn, 0, 0, g{fi), 0, /i), (4.72) 

for K 0. Here ci is a dimension of correlation function and 

r'M=exp[- r^7(^7(/),0)]. (4.73) 
J fj,' fi 
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A A Note on the Combinatoric of Planar Graphs of 
a Matrix Model 



In this Appendix we prove the following statement. The planar correlation functions 
without vacuum insertions in all orders of perturbation theory arc represented by the 
sum of all topologically non-equivalent graphs without any combinatiric factors. 

The proof of this statement follows from the following considerations. Let us consider 
an arbitrary planar graph with n external lines and with m four-point vertices. We will 
consider the external lines of graph corresponding to global invariant Green functions as 
the lines corresponding to an n-point generalized vertex. We suppose that ki from m 
four-point vertices are connected with the generalized vertex at least by one line (see Fig. 
16, where all double lines are drawn by single lines for simplicity). 

The combinatoric factor associated with all possible contractions of ki vertices with 
the n-point generalized vertex is equal to 

(m — ki}\ 

Indeed, to produce topologically equivalent graphs ki vertices can be chosen from m 
vertices by ^jn^[yj,^\ ways, ki vertices can be permutated by ki\ ways and each vertex can 
be attached to the generalized vertex by 4 ways. 

After connecting the generalized vertex by lines with ki vertices we get the graph 
which can be considered as some set of generalized vertices which are connected with 
m — ki remaining four-point vertices. For example, on Fig. 16 there are at least two new 
generalized vertices. They are formed by lines outcoming from 4- vertices labeled by 1, 2, 
3, 6, ... and 4, 5. Let us supposed that there are r new generalized vertices and that k2 
from m — ki remaining four-point vertices are connected with new generalized vertices at 
least by one line so that the first generalized vertex is connected with k2'^ 4-vertices, the 
second generalized vertex is connected with 4-vertices and so on. There are 

(m-fci)! fca! ,(1), ,>),_ (m-fci)! . 

{m-k,-k2)\k2\ ki'^\...k^;h ' ■■■■ ' ■ {m-ki-k2)\ 

possibilities to do this. Again we get the graph that can be considered as some set of 
generalized vertices which are connected in some way with m—ki—k2 remaining four-point 
vertices. 

Analyzing the construction of a full graph by steps that consist in attaching some 
vertices to generalized vertices and assuming that it is need to do I steps to get the graph 
one gets that the total combinatoric factor will be equal to 

(m — ki)\ {m — ki — A;2)! {m — ki — k2 — ... — ki)\ 

This factor is just canceled by the factor l/(m!4"*) appearing from the expansion of the 
exponent of the interaction Lagrangian, that gives the prove of the statement. 
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Figure 1: Derivation of the Schwinger-Dyson equations for free Green's functions in the 
Boltzmannian Fock space. 
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Figure 2: Derivation of the Schwinger-Dyson equations for free Green's functions in the 
Boltzmannian Fock space (the 2A;-th marked operator). 
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Figure 3: Derivation of the Schwinger-Dyson equations for free Green's functions in the 
Boltzmannian Fock space (the 2k + 1-th marked operator). 
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Figure 4: Derivation of the Schwinger-Dyson equations in the Boltzmannian Fock space 
for the interaction 0*". 
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Figure 5: Planar graphs in the matrix theory with the quartic interaction up to order g^. 
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Figure 7: Half-planar graphs in the Boltzmann theory with the quartic interaction in the 
order g^. 
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Figure 8: The dependences of the two-point Green's functions Gn and Gu on the couphng 
constant g for different values of 




Figure 9: Graphical representation of equation (4.21). 
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Figure 10: Graphical representation of equation (4.20). 




Figure 11: Graphical representation of equation (4.23). 
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Figure 12: Graphical representation of equations (4.29) and (4.30) 
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Figure 13: a) Example of divergent parts of Il-type of the two-point half-planar graph, 
b), c) Examples of half-planar divergent subgraphs. 
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Figure 14: Two possibilities making a contraction of Il-type subgraph to a point. 
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Figure 15: Graph for 4-point correlation function in the Boltzman field theory. Some 
fines of grapfi a) are marked by crosses. Insertions of SM'^cf)'^ and {Z^ — 1)((9^0)^ are 
admissible only on the marked lines. 



47 



Xi Xi X-i X4 



Figure 16: Construction of planar graph, ki of m four-point vertices are connected with 
the generahzed vertex at least by one line. Two new generalized vertices are formed by 
hnes outcoming from 4- vertices labeled by 1, 2, 3, 6, ... and 4, 5. 
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